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ABSTRACT 

A new technique of Fourier analysis on a cut sky (Gorski, 1994) has been 
applied to the two year COBE DMR sky maps. The Bayesian power spectrum 
estimation results are consistent with the Harrison-Zel'dovich n = 1 model. 
The maximum likelihood estimates of the parameters of the power spectrum 
of primordial perturbations are n = 1.22 (1.02) and Q rm s-PS = 17 (20) fiK 
including (excluding) the quadrupole. The marginal likelihood function on n 
renders n = 1.10 ± 0.32 (0.87 ± 0.36). 

Subject headings: cosmic microwave background — cosmology: observations 
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1. INTRODUCTION 

The determination of the power spectrum of primordial inhomogeneities and their 
consistency with the predictions of inflation are critical issues in contemporary cosmology. 
Standard inflationary scenarios imply a power law spectrum P(k) oc k n } with n ~ f on 
the scales probed in the COBE DMR sky maps. Previous attempts to determine the 
primordial power spectrum from the DMR anisotropy data (Smoot et al. 1992, Wright 
et al. 1992, Adams et al. 1992, Scaramella & Vittorio 1993, Seljak & Bertschinger 1993, 
Bennett, et al. 1994, Smoot et al. 1994, Wright, et al. 1994a, Bond 1994) found both 
relatively steep spectra and sensitivity to the inclusion or exclusion of the quadrupole in 
their results. These methods have all employed approximate statistical techniques and have 
relied on Monte Carlo techniques to assess and/or calibrate the final results. 

In this Letter, the COBE DMR two year 53 and 90 GHz data are analysed to 
determine the primordial power spectrum using the method described in Gorski (1994). 
Specifically, the sky maps are Fourier decomposed in the basis of orthonormal functions 
on the Galaxy cut sky to yield a set of harmonic mode coefficients. These are then used 
in a maximum likelihood analysis to infer the parameters of the theoretical anisotropy 
models. The merits of the present method are: 1) harmonic mode coupling is explicitly 
accounted for by construction of the orthonormal functions on the Galaxy-cut sky, 2) since 
the harmonic modes have a Gaussian probability distribution (Gorski, 1994) an exact 
likelihood function for the model parameters can be employed, 3) the technique permits a 
simultaneous analysis of different frequency maps taking full advantage of both the auto- 
and cross-correlation information in the data. 

In this analysis we Fourier decompose the two year DMR sky maps over the spectral 
range up to £ = 30, where both the DMR beam response has fallen to ~ 0.2 (Wright et 
al. 1994b) and the multipole amplitude is noise dominated. The monopole and dipole 
components, which are physically irrelevant for the power spectrum estimation, are 
algebraically excluded (Gorski 1994). In addition, no attempt is made to formally model 
and subtract the high latitude galactic emission, which is predominantly quadrupolar in 
nature. Therefore, in what follows, the power spectrum parameters are derived for two 
cases: one in which data spanning the multipole range £ £ [2, 30] is used, and the other in 
which the quadrupole (l = 2) mode is excluded. 

Hereafter, bold, upper case letters denote matrices; bold, lower case letters denote 
vectors; and p is a pixel label. Specific equations from Gorski (1994) will be referred to by 
G immediately followed by the equation number. 
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2. DATA 

Each individual DMR sky map is a collection of pixelized measurements of temperature 
perturbations, A(p), comprising cosmological, receiver noise, galactic, and systematic 
signals. The instrument noise is characterized by the number of observations per pixel, and 
an rms noise per observation, <t (derived from Table 1 of Bennett et al. 1994). In order 
to limit the effects of galactic emission on the power spectrum determination, a \b\ = 20° 
Galaxy cut is imposed on the maps. Systematic effects have been strongly limited (Bennett 
et al. 1994) and are neglected in the present analysis, as are the weak high latitude galactic 
signals. The 53 and 90 GHz two year cut-sky maps are analyzed in this work. We form 
weighted average maps at each frequency by combining the individual A and B channels as 
follows 

A 53 = 0.578 A 53 a + 0.422 A 53B , 

A 90 = 0.366 A 9(L4 + 0.634 A 90B , (1) 

where all temperatures are in thermodynamic units. These weights were chosen to minimize 
the average noise variance per pixel in the combined maps. Fourier analysis of the cut-sky 
maps is performed in the 961-dimensional linear space, spanned by the orthonormal 
functions if) constructed in Gorski (1994). The harmonic coefficients of a map in the if) 
basis are given by (see G8 and G3) 

Mp)\ P e{cut sky} = c' • V> (p), where c = (A tf) ) {cut sky} . (2) 
The data vector for joint analysis of the 53 and 90 GHz data was formed according to 

C T _ / X X \ 
*-53ffi90 — l c 53' '-90/' • 

3. LIKELIHOOD ANALYSIS 

The covariance matrices, C = Ccmb + Cjv, of the harmonic amplitudes, c = Ccmb + Cjv, 
specify the Gaussian probability distribution (G15) of the receiver noise contaminated 
theoretical CMB anisotropy signal. The power spectrum dependent matrix Ccmb can be 
conveniently parameterized by two quantities, Q rm s-PS and n (see G10, and Gil). 

The receiver noise is modeled as a spatially uncorrelated Gaussian random process 
with a frequency dependent variance per pixel a 2 (p) = ( A 2 N (p)) {nmse ensemble y = a%/N(p). 
Eq. G13 was used to obtain the Fourier space noise correlation matrices Cjv- The resulting 
rms noise amplitude per harmonic mode is 4.8//K for the 53 GHz data, and 7.8//K for the 
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90 GHz data. These values are important for an a priori assessment of the sensitivity of 
DMR to the anisotropy predicted by a particular cosmology. 

Having established an appropriate mathematical representation of the data, specifically 
the Fourier amplitude vector 053990, the anisotropy correlation matrix Ccmb and the noise 
correlation matrices Cjv53 and Cjv9o, we can perform a Bayesian likelihood analysis to 
determine the power spectrum parameters Q rm s-PS and n. The necessary ingredients for 
a proper Bayesian statistical analysis (e.g. Berger 1980) of values of Q rm s-PS and n to 
be inferred from the data are the likelihood function, P(c 53(S90 \Q rms _ps } n), and the prior 
distribution Tr(Q rms _ps } n). These then define the posterior distribution 

P(Qrms-PS, ™|c 5 3ffi9o) ^ ^(Qrms-PS, ™) X P(c 

53090 I rms — ps,n). (3) 

We evaluate the likelihood function for Q rm s-PS and n (see eqs. G15 and G16) given data 
either including or excluding the quadrupole anisotropy (which involves correlation matrices 
of dimension 1914, or 1904, respectively). 

Figure 1 shows two views of the likelihood function evaluated with the quadrupole 
anisotropy included; the no quadrupole case yields qualitatively similar contours. This 
likelihood function is confined to a very narrow and steep ridge in the parameter space 
conveniently defined by the multipole amplitude a 9 = 8 fiK. Clearly, any prior Tr(Q rms _ps } n) 
that is slowly varying near the peak of the likelihood will not significantly affect the inference 
on Qrms-PS and n. We have tested the sensitivity of the derived maximum likelihood 
parameter values and credible intervals to some specific assumptions about tt. The following 
prior distributions were considered: (1) tt = const, (2) 7r oc 1/(1 + n), (3) 7r oc (1 + n), (4) 

7T OC l/Q r ms-PS, (5) 7T OC Q rm s-PS, (6) 7T OC 1 / ( 1 + n) /Qrms-PS , and (7) 7T OC (1 + Tl) Q rm s-PS- 

As expected the results are robust to ~ 1 //K in Q rm s-PS^ and ~ 0.1 in n. Hence, in the 
following we adopt a constant prior distribution. 

Figure 2 shows the 68%, 95%, and 99.7% contours for the posterior distributions of 
Qrms-PS and n evaluated with and without the quadrupole anisotropy. Exclusion of the 
quadrupole from the analysis results in considerable elongation of the likelihood contours 
towards low n and high Q rm s-PS values. This is the region of the parameter space seriously 
constrained by the inclusion of the relatively low observed quadrupole (Bennett et al. 1994). 
Marginalization by integration of the likelihood function over Q rm s-PS renders a likelihood 
density for the spectral index n shown in Fig. 3. The conditional likelihood density on 
Qrms-PS for n = 1 (the Harrison-Zel'dovich spectrum) is shown in Fig. 4. Relevant numbers 
are presented in Table 1. 



4. DISCUSSION 
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Miscellaneous issues pertaining to the proper assessment of our results include (1) noise 
model uncertainties, (2) the question of biases in parameter inference, and (3) comparison 
to other attempts at determination of the power spectrum. 

The method of power spectrum determination presented in this Letter requires an 
adequate description of the statistical properties of noise in the sky maps. In order to 
test the correctness of construction of the Gaussian noise ensemble, we have evaluated the 
Fourier components of the cut-sky two year difference maps, Q-a-b^ at both frequencies, and 
the a priori difference map noise covariance matrices (see eq. GI3), Cjv(a-B)- The reduced 
X 2 , defined as c\_ B ■ Cn(a-b) • ca-b/N where N = 957, is 0.93 at 53 GHz and 1.06 at 90 
GHz. Therefore, the a priori noise model appears sufficiently accurate when constructed 
using the flight rms noise values. Nevertheless, the noise rms per observation, <t , is only 
known to an accuracy of ~ 1%. The possibility of a mismatch between the assumed and the 
actual noise levels can perturb the spectral index and amplitude estimation. Specifically, 
an underestimate of the noise leaves some small scale power which the likelihood method 
translates into a steeper spectrum than required by the true sky anisotropy, or vice versa. 
If cr is adjusted by ~ 1%, the maximum likelihood values for Q rm s-PS and n are shifted 
by ~ 0.5 fiK and ~ 0.04, respectively, insignificant amounts when compared to the overall 
uncertainty in the inferred spectrum. 

In principle, maximum likelihood estimates of non-linear model parameters are only 
asymptotically unbiased. In order to test whether any biases exist in this analysis, Monte 
Carlo simulations of noisy random skies drawn from an ensemble of fixed Q rm s-PS and n 
and with the noise properties of the two year 53 GHz maps were analysed. No statistically 
significant bias in the simulated sample-averaged estimates of the parameters Q rm s-PS and 
n was detected when the quadrupole was either included or excluded. A sub-ensemble of 
simulations with suppressed quadrupole amplitudes was also considered: in this case, a 
general bias was observed to lower Q rm s-PS and higher n. Given that the observed sky 
quadrupole is low with respect to our inferred ensemble averaged value of Q rm s-PS = 20 //K, 
this affords a plausible explanation for the small Sn ~ 0.2 difference between spectral slopes 
inferred from the sky maps with or without the quadrupole (see Table I). 

By choosing £ m ax = 30 we have analysed the COBE DMR maps over a spectral range 
larger than previously implemented. Other methods are inherently insensitive to the high-£ 
information content of the data either by construction (e.g. Wright et al. 1994, where 
^max = 19 was adopted), or as a consequence of additional binning or smoothing of the 
data (e.g. Bennett et al. 1994, Smoot et al. 1994). Specifically, we have tested our method 
over the restricted range £ £ [3, 19] as a direct comparison with Wright et al. (1994). The 
maximum likelihood estimates in this case are Q rm s-PS = 17.8 //K, and n = 1.22. The 
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steepening of the resulting spectrum, in accord with Wright et al. (1994), is driven by the 
relatively high amplitude of several modes within £ £ [14, 19] in the 53 GHz sum map. In 
addition, we find that the power in the corresponding bin of the 53 GHz (A-B) difference 
map falls below the a priori noise model, thus Wright et al. determine a still higher power 
in their quadratic power spectrum estimator. With the inclusion of £ > 20 modes the 
effect of this excess is more effectively constrained and attributed to noise. A more detailed 
comparison of the application of different methods of power spectrum inference is deferred 
to a future paper. 

In summary, the new, statistically unbiased technique of power spectrum determination 
described in Gorski (1994) has been implemented for the COBE DMR 53 and 90 GHz 
two year sky maps. The results are completely consistent with the Harrison-Zel'dovich 
n = 1 spectrum favored by inflationary models. If this is indeed the case, then the two 
year COBE DMR data provide a ~ 12a significant detection, Q rm s-PS = 20 //K, for the 
amplitude of the spectrum of primordial inhomogeneity. 

We gratefully acknowledge the efforts of those contributing to the COBE DMR. 
COBE is supported by the Office of Space Sciences of NASA Headquarters. 
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Table 1: Derived parameters from the likelihood analysis 





including Quadrupole 


excluding Quadrupole 


Qrms-PS (/^K) 


n 


Qrms-PS (/^K) 


n 


Maximum likelihood values 


17.0 


1.22 


20.0 


1.02 


68% credible interval 


[12.2, 24.6] 


[0.70, 1.68] 


[13.5, 30.5] 


[0.43, 1.55] 


95% credible interval 


[10.0, 31.8] 


[0.35, 1.95] 


[10.6, 41.5] 


[0.02, 1.85] 


Marginal likelihood on n 


1.10 ± 0.32 


0.87 ± 0.36 


Conditional likelihood on Q^Zl-PS 


19.9 ± 1.6 


20.4 ± 1.7 
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Fig. 1. — Two views of the likelihood function, P(n } Q rms _ps) } derived in a simultaneous 
analysis of the 53 and 90 GHz COBE-DMK two year data including harmonic amplitudes 
from the range £ £ [2, 30]. 
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Fig. 2. — Contour plots of the likelihood function P(n } Q rms _ps). 68%, 95% and 99.7% 
levels and the maximum values are shown. Solid lines and ' + ' represent the results of the 
analysis including the quadrupole. Dotted lines and '*' show the no quadrupole case. 
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Fig. 3. — Marginal likelihood densities, p(n) oc / dQ rms _ps P(n } Q rms _ps) } for the spectral 
index n (arbitrary normalization). The solid curve includes the quadrupole in the analysis, 
the dotted curve excludes the quadrupole. 
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Fig. 4. — The n = 1-conditional likelihood density (arbitrary normalization) for Q rms _ps. 
The solid curve includes the quadrupole, the dotted curve is the no quadrupole case. 



